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Qualitative properties of coexistence and semi-trivial limit 
profiles of nonautonomous nonlinear parabolic Diriclilet 

systems 
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Abstract 

We study the symmetry properties of limit profiles of nonautonomous nonlinear parabolic 
systems with Dirichlet boundary conditions in radial bounded domains. In the case of 
competitive systems, we show that if the initial profiles satisfy a reflectional inequality with 
respect to a hyperplane, then all limit profiles are foliated Schwarz symmetric with respect 
to antipodal points. One of the main ingredients in the proof is a new parabolic version 
of Serrin’s boundary point lemma. Results on radial symmetry of semi-trivial profiles are 
discussed also for noncompetitive systems. 
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1 Introduction 


The main motivation for this paper is the study of the asymptotic shape of positive global 
solutions of the following nonautonomous nonlinear parabolic problem. 

{ui)t - Aui = fi(t, \x\,Ui) - ai{\x\,t)uiU2, x C B, t > 0, 

Ui{x,t) = Q, xCdB, t > 0, ( 1 . 1 ) 

Ui(x,0) = uoj(x), xeB 


for i = 1,2. Here and in the remainder of the paper, B denotes a ball or an annulus in 
with N > 2. The nonlinearity / and the coefficients cti satisfy some smoothness assumptions. 
As an important example of a problem to which our results apply we have the competitive 
Lotka-Volterra system, that is, 


(ui)t - Aui = ai{t)ui - bi(t)ul - ai{t)uiU2, 
{U2)t - Au 2 = a2{t)u2 - b2{t)ul - a2{t)uiU2, 
Ui = U2 = 0 
Ui{x, 0) = uo,i{x) > 0 


X G B, t > 0, 

X € B, t > 0, 
on dB X (0, oo), 
for X € B, i = 1, 2, 


where for some /3 £ (0,1) and i = 1,2, we have that 

Oi, bi, G oo)) n L°°((0, oo)) and 


inf ai{t) > 0. 


( 1 . 2 ) 


(1.3) 
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For a discussion on the history and the relevance of system (EH) we refer to HHl] and the 
references therein. 

To describe the asymptotic shape of a solution u = (mi,U 2 ) of (11.21) we use the omega limit 
set, that is, 

uj{u) =u){ui,U 2 ) ■■= {{zi,Z 2 ) & C{B) X C{B) : (1.4) 

lim \\ui{-,tn) - + ||u 2 (-,t„) - Zi||L“(B) = 0 for some -)► oo}. 

n—^oo 

This set is nonempty, compact, and connected for uniformly bounded solutions of (EH) , see 
Lemma |3.II below. Elements of uj{u) are called limit profiles. If z = (zi,Z 2 ) G u}(u) is such that 
zi ^ 0 and Z 2 ^ 0 in B, then we call z a coexistence limit profile. If Zi = 0 and Z 2 ^ 0 or if 
Z 2 = 0 and zi ^ 0, then we call z a semi-trivial limit profile. In this paper we present results on 
the symmetry of both kinds of limit profiles. The proofs rely on a variant of the moving plane 
method —often called rotating plane method —and on stability arguments. 

Let us introduce some notation that will be used throughout the paper. Let = {a; g 

: |x| = 1} be the unit sphere in N > 2. For a vector e € we consider the 

hyperplane H{e) := {x € R^ : a: • e = 0} and the half domain B{e) := {x G B : x ■ e > 0}. We 
write also ae ■ B ^ B to denote reflection with respect to H{e), i.e. a-e{x) := x — 2{x ■ e)e for 
each X G B. Following m, we say that a function u G C{B) is foliated Schwarz symmetric with 
respect to some unit vector p G if u is axially symmetric with respect to the axis Mp and 

nonincreasing in the polar angle 0 := arccos( j||- • p) G [0, tt]. 

Our main result for system E3 is the following. 

Theorem 1.1. Suppose that (11.31) holds and let u = (mi,U 2 ) be a classical solution of (11.21) such 
that Ui G L°°{B x (0, oo)) for i = 1,2. Moreover, assume that 

I for some e G with 

(hO) uo,i > ito,i O (Te, uo, 2 < Mo,2 o (Te in B{e) < „ 

[ Mo,i ^ Mo.i o cTg for 1 = 1,2. 

Then there is some p G such that all elements (zi,Z 2 ) G lo(u) satisfy that zi is foliated 

Schwarz symmetric with respect to p and Z 2 is foliated Schwarz symmetric with respect to —p. 

This theorem is an direct consequence of the next more general result. In the following 
(^a,a/ 2 (Q) for a > 0 denotes the standard (parabolic) Holder space on Q C see (13.11) . Fix 

Ib ■■= {|a;| : x G B}. 

Theorem 1.2. Let ui,U 2 G C^'^(B x (0, oo)) (~) C{B x [0, oo)) (~) L°°[B x (0, oo)) be nonnegative 
functions such that u = {ui,U 2 ) is a solution of system (II. 1|) where the following holds. 

(hi) For i = 1,2, the nonlinearity fi : [0, oo) x x [0, oo) — > R, {t, r, v) i—>• fift, r, v) is continu¬ 
ously differentiable in v. Further, the functions fi and dyfi are Holder continuous in t and 
r for all v G [0,oo), and locally Lipschitz continuous in v uniformly with respect to t and 
r. In other words, there is 7 > 0 such that, for every h G {fi, dyfi : i = 1, 2}, 

h(-, ■, v) G X (0, 00 )) for all v G [0, 00 ) 


and 


\h{t,r,v) - h{t,r,v)\ 

I -1 

rG/B,i>0, P ~ '^1 

v,v^K,v^v 


< OO 


for any compact subset K C [0,oo). Moreover fi(t,r,0) = 0 for all r G Ib, t > 0, and 
iG {1,2}. 
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(h2) There are positive constants a*, a*, and fi sueh that ai € x (0,oo)) and a, < 

ctiir, t) < a* for all r € Ib, t > 0, and i £ { 1 , 2 }. 

(hS) There is e G for which uo,i ^ uo,i o o-g, uq ,2 ^ uo ,2 o cce, and uq^i > uop o Cg, uo,i < 

Wo ,2 o CTg in B(e). 

Then there is some p £ such that all elements {zi,Z 2 ) £ w(u) satisfy that zi is foliated 

Schwarz symmetric with respect to p, and Z 2 is foliated Schwarz symmetric with respect to —p. 

Theorem ll.2l is a Dirichlet analog of [T^ Theorem 1.3], where the asymptotic shape of solutions 
of the Neumann version of system m is studied, and also foliated Schwarz symmetry with 
respect to antipodal points is proved. In both theorems assumption (h3) is essential to perform 
a rotating plane method and the symmetry result can not hold in general without it; see m 
and m for a discussion in this regard. 

One of the main ingredients in the proofs in m is an extension of the solution using a 
reflection with respect to the boundary of the domain. This extension is only possible under 
Neumann boundary conditions, and in fact, the main perturbation tool |161 Lemma 3.1] can not 
hold under Dirichlet boundary conditions. 

In this paper, to prove Theorem ll.2l we use a new parabolic version of Serrin’s boundary point 
lemma m Lemma 1] that provides quantitative information on the second derivatives at corner 
points. This information is used to perform a perturbation argument within a rotating plane 
method scheme. A parabolic version of Serrin’s boundary point lemma was obtained in [U Lemma 
2.1], where the result from Serrin nzi Theorem 1] is extended to the parabolic setting; however, 
this version is insufficient for our purposes, since it is crucial to obtain estimates independent of 
a specific solution, see Lemmas 12.31 and 12.41 below. 

In this setting, to perform the perturbation argument one needs some equicontinuity of the 
solution up to second derivatives. This is the reason why the right hand side in and in (ED 
is required to satisfy some Holder regularity (this assumption is not needed in [16) . but it is a 
common hypothesis in the study of nonautonomous Lotka-Volterra models, see [U Section 5.5]). 
With this hypothesis one can guarantee using some standard regularity techniques as in m 
and [2] that the solutions ui and U 2 belong to for some 7 £ (0,1) (see Lemma 

13.11 belowl. 

As already noted in [16] , the possibility of the existence of semi-trivial limit profiles presents a 
complication for the perturbation argument, which depends on lower bounds for some norms 
of solutions to a linearized problem (see Proposition 12.71 belowl. In [TB] this complication was 
circumvented using a normalization argument relying again on estimates which only hold for 
the Neumann problem. Here we adapt this argument to the Dirichlet case using some sharp 
estimates due to Hiiska, Polacik, and Safonov for positive solutions of homogeneous linear 
problems. 

For a class of parabolic cooperative systems, Fdldes and Polacik [5] proved in particular that, 
if the domain is a ball and the nonlinearity satisfies some monotonicity assumptions, positive 
solutions are asymptotically radially symmetric and radially decreasing. Their proofs rely on 
maximum principles for small domains, which is a different approach from ours and can not be 
applied here directly (see also Remark 14.41 in this regard). See also [15] for a version of Theorem 
O for more general scalar equations using maximum principles for small domains. 

As explained in II US], systems of the type (ED are commonly used to model population 
dynamics. In this setting, the existence of a semi-trivial limit profile can be interpreted as 
the asymptotic extinction of one of the species. In such a situation, one may guess that the 
population density of the remaining species, in the (asymptotic) absence of a competing (or 
symbiotic) species, is likely to become asymptotically radially symmetric in B. We show in |14j 
that this is indeed the case if additional assumptions are made. 
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Let Ai > 0 denote the first Dirichlet eigenvalue of the Laplacian in B. 

Theorem 1.3 (Theorem 5.2 in [2]). Let ui,U 2 G x (0, oo)) fl C{B x [0, oo)) fl L°°{B x 

(0,oo)) be nonnegative functions such that u = (ui,U 2 ) is a classical solution of 

{ui)t - Aui = fi{ui) - ai{x, t)uiU 2 in B X (0, oo), 

Ui =0 on dB X (0, oo), 

where m(-,0) = G C{B) is not identically zero, at G L°°{B x (0,oo)), fi G C'^(R) is strictly 
concave in [0, oo), fi(0) = 0, and lim fi{s) = —oo for i = 1,2. If /((O) > Ai and /^(O) > Ai, 

s—^oo 

then all the semi-trivial limit profiles of u are radially symmetric. In particular, if there is 
(z,0) G w(m), then z is radially symmetric and it is the unique positive solution of —Az = fi{z) 
in B. The analogous claim holds if there is (0, z) G uj(u) with z ^ 0. 

The proof relies on well-known techniques for nonlinear PDE’s with concave nonlinearities 
as used in [7] (see also [3 Chapters 9,10]). The hypothesis of Theorem 11.31 might look a bit 
restrictive, but they are satisfied by many interesting models; for example, the Lotka-Volterra 
models for two species. We state this next. 

Corollary 1.4. Let ui,U 2 G x (0, oo)) fl C{B x [0, oo)) fl L°°{B x (0, oo)) be nonnegative 

functions such that u = (mi,M 2 ) is a classical solution of 

{ui)t — Aui = QiUi — biU^ — ai{x, t)uiU 2 in B x (0, oo), 

Ui = 0 on dB X (0,oo), 

where u{-, 0) = itoy G C{B) is not identically zero, Oi > Ai, bi > 0, and ai G L°°{B x (0, oo)) for 
i = 1,2. If {z, 0) G uj{u) then z G C{B) fl Hq{B) is a positive radially symmetric function and it 
is the unique solution of —Az = aiz — biz^ in B. The analogous claim holds for (0,z) G uj{u). 

The Lotka-Volterra model for competing species assumes that the coefficients ai are positive 
for i = 1,2, but this assumption is not needed for the claim in Corollarv ll.41 Indeed, Theorem 1 1.3 1 
makes no assumption on the sign of a\ and 02 , and therefore the result also applies to semitrivial 
limit profiles of cooperative, competitive, and mixed type systems, which includes, for instance, 
the predator-prey model. Also note that the coefficients ai may depend on x and not only on 
\x\ as in Theorem 11.21 Te keep this paper short, we do not give here the proof of Theorem 11.31 
We refer to m for the proof and other related results. 

We now give a brief outline of this paper. Section [2| is devoted to the parabolic version 
of Serrin’s boundary point lemma for supersolutions of scalar equations and its application to 
prove a perturbation lemma for solutions to weakly coupled linear systems. In Subsection 12.II we 
recollect some estimates from [3, which will be used for the normalization argument. In Section 
13 we explain how the equicontinuity up to second derivatives of any solution is achieved and 
finally Section contains the proof of Theorem 11.21 

To close this introduction, let us mention that, similarly as in |16j . the proof of Theorem 1 1.2 1 
can be adjusted to deal with other kinds of systems. For example, cubic systems, a cooperative 
version of (lEa), or irreducible cooperative systems of n—equations can be considered. In this 
last type of systems, no normalization procedure is necessary and therefore also sign-changing 
solutions are allowed. For the linearization of these systems one can use the Hadamard formulas 
similarly as in [B]. Then it is relatively easy to adjust the rotating plane method presented here 
to these problems. Since the statements and the proofs of these results are similar to those 
presented in m and in m for the Neumann case, we do not give any further detail in this 
regard. 
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2 A parabolic version of Serrin’s boundary point lemma 

First we fix some notation. Let B he a. ball or an annulus in centered at zero and fix 
O<A 1 <A 2<00 such that 

^ _ f{a; € : Ai < |a:| < A 2 }, if Ai > 0, 

|{a; € : |a:| < A 2 }, if Ai = 0. 

We denote by Br{y) C R^ the ball of radius r centered at y G R^. For a subset A C B and 
(5 > 0 we define 


[A]^ := {x G B : dist(a;, A) < 5}. (2.2) 

For X C R, e £ and a function v : B x I ^Rwe define 

v^:BxI^R by v^{x,t)'■= v{x,t) — v{ae{x),t)- 


Fix I := [0,1], ei = (1,0,..., 0) G R'^, and let r; : B x J —>• [0, 00 ) be a function such that 
satisfies 

>0 in B(ei) x /, 

=0 on dB{ei) x I, (2.3) 

>0 in B(ei) x I, 

and the following holds. 

{Ha,Po) There are a G (0,1) and /3o > 0 such that 


\v\2+a-,B^I + ||c||L~(B(e^)x7) < /3o, 
where | • | 2 +a;Sx/ denotes the parabolic Holder norm, see (13.11) . 


[Hk) There is fc > 0 such that \\v‘^^\\L'=°{B{e)x{^A)) — 

Remark 2.1. Let v G C^’^(B x I) be a function satisfying and {Hk). If = 0 on 

dB{e) X I, then there is (5 £ (0, ); depending only on a, /3o, k, and B such that \ < k 

in [dB{e )\5 x I. Therefore we have that |lLo=((B(e)\[as(e)],)x(A|)) > 

We will use the following version of the parabolic Harnack inequality for general domains as 
given in m- 

Lemma 2.2 (Particular case of Lemma 3.5 in m)- Let /3o > 0, £ > 0, J > 0, 0 > 0, 0 < r < 
Ti < T2 < T3 < T4 < 00, be some given constants. There are positive constants k, ki, and p 
determined only by B, /?o, S, e, T2 — T1, T3— T2, T3 — T4, and 6 with the following property. If D, U 
are subdomains in B with D CC U, dist{D,dU) > S, \D\ > e, and v G loci^ ^ bl 

C{U X [ t , r4]) is such that v > 0 in U x ( t , T4) and satisfies Vt — Av — cv > 0 in U x (r,T4), 
where ti — 20 < t < ti — 6 and ||c||itxj(f/x(r,T4)) < Po. Then 


inf V > K 

Dx{t3,T4,) 


1 

\D X (ri,T2)| 


f yP 

J Dx{ti,T2) 
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Indeed, this follows directly by m Lemma 3.5] taking v > 0 and g > 0. 

The following lemma focuses on points of the boundary dB[ei) which are not corner points. 

Lemma 2.3. Let v € x /) be a function such that satisfies (12.311 and such that 

assumptions {Ha^pg) and {Hk) hold. Then given 6 € (0, ), there are positive constants e 

and g depending only on <5, a, /3o, k, and B such that 

v‘^^{x,t) > gxi for all X G B{ei)\[dB]s, (2.4) 

and 

> e in {dB{ei)\[dB r\ H{ei)]s) X [^ , 1], (2.5) 

Here v is the inwards unit normal vector field on dBici). 

Proof. By Remark I^T] there is £ (0,5) depending only on a, /?o, k, S, and B, such that 

lk'''llL~((B(e)\[aB(e)],i)x(i.|)) > (2-6) 

Let r < min{5i, j} be such that for any t* £ [|, 1] and x* £ dB{ei) with dist(a:*, iL(ei)(~l9i3) > 5i 
we have that x* +rh'{x*) =: y £ B{ei) and dBr{y,t*) (H {dB{ei) x [0,t*]) = {{x*,t*)}. Fix such 
a pair {x*,t*) and define the sets 

D:=Br{yX)r)Br(x*,t*)niB{ei) x [0,^]), 

Fi := dBriy, t*) n Br{x*,t*) n (B{^ x [0, t*]), 

Fa := Briy,t*)ndBr{x*,t*) n (B{^ x [0,^]), 
and the function 2; : D —[0,1] given by 

z{x, t) := _ e-T'’-')e-^'>‘ with 7 = 

An easy calculation shows that 


zt — Az — cz 

— g-7(b-y| +(i-t ) )-Pot2j(^_2j\x — y\^ + N — (t — t*)) + (—/3o — c)z 
< g-7(b-yP+(i-t*)")-/3ot^('_^^2_^2(iV + l)) =0 inL>. (2.7) 

Let 0 < d < min{dist(F2,9i3(ei) x [|,1]) , 5i} and K := {x £ B{ei) : dist(a;, 9i3(ei)) > d}. 
Note that r2 C AT x [|, 1] and K only depends on a, /3o, k, 6, and B. Then by {Ha,pg), (12.6|) . and 
Lemma there is gi = g,i{a,l3o,S,k,B) > 0 such that > gi in K x [|, 1]. In particular, 
we have that v^^ > /ii in r2. Let w : D ^ M. he given by w := — giz. Since z = 0 in Fi 

and > 0 in B{ei) we have that i/; > 0 in Fi and, since z < I in F2, it follows that w > 0 
on Fi U r2. Further, by (12.31) and (12.71) we get that Wt — Aw — cw >0 in D. Thus the parabolic 
maximum principle implies that w > 0 in D. Since w{x*, t*) = 0 we have that ^ > 0 and 

thus 

(x*,t*) > gi^{x*,t*) > 2e~'^^^~^°gi-fr -. e > 0. 
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Since 5 > (5i, this proves (12.51) . Moreover, since > e in {H{ei)\[dB]si) x there is, by 

{Ha,po), some S 2 = 62 {a,l 3 o,S,k,B) G (0, ^i) such that 

^ > I in {[H{e,)\[dB]s,]s,) x [^, 1]. (2.8) 


On the other hand, i?(ei)\[9i3(ei)]5j C B{ei)\[dB{ei)]s 2 since S 2 < ^i- Thus, by (Ha^po), 
(ESI), and Lemma l2.21 there is /r2 > 0 depending only on a, /3o, <5, fc, and B such that > fi2 in 
(i3(ei)\[9B(ei)]52) X [|, 1]. This, together with (12.81) easily implies (12.4|) for some /i > 0 depending 
only on a, /3o, k, and B. □ 

We now turn our attention to the comer points on the boundary dB{ei). 

Lemma 2.4. Let v € x /) be a function such that satisfies (12.31) and sueh that 

assumptions and (Hk) hold. Then there is e > 0 depending only on a, j3o, k, and B such 

that 


1 ) 

ds"^ 


> e 


and 


(x, 1) 

dP 


< —e 


(2.9) 


for all X G dB O H{e-i), where s = -^( j ^ + ei) G s = + ei) G and v is the 

inwards unit normal vector field on dB. 

Proof. Let x* G dB O H{ei) and 0 < r < min{i, }, see (12.11) . Then x* + rv{x*) =: y G 

H{ei) n B and dBr{y, 1) (~l {dB x I) = {(x*, 1)}. Define 


U := Briy, 1) n Sr (x*, 1) O (B(ei) x J), 

Ai := dBr{y, 1) (~l Br{x*^ 1) 0 {B{ei) x /), 

A 2 := Br{y, 1) n dBr{x\ 1) O (B(^ x J), 

A3 := Br{y,l) n Br(x*,l) n (iJ(ei) x /), 

and the function (p ■. B x I ^ [—A2, A 2 ] given by 

p{x,t) ■■= - e-®"')e-'^'>* with 9 = . 

A direct calculation shows that 


Pt - - cp 

= e-®ili^’‘i-(^’i)l'i-'^“‘26(xi(-26(|x - y\^ + N + 2 - t) - {^o + c)p 
< g-e(lfo,t)-(y.i)l")-/3ot25/^^(_26»(^)2 _p^ + 3) = 0 in D. (2.10) 

Since dist(A 2 , dB x /) > 0, there is = Si{B) > 0 such that A 2 C (B(ei)\[ 9 i 3 ] 5 J x [ |, 1 ]. 
Then, by Lemma 12.31 with S = i5i, there is /i > 0 depending only on a, /3o, k, and B such 
that x®i(x,t) > ^Xi for x G {B{ei)\[dB] 5 ^) and t G [|, 1]. In particular v^^{x,t) > yxi for 
(x,t) G A 2 . Define the function ^ B x / —>• R by 'if{x,t) ■= v^^{x,t) — Then ip > 0 

on A 2 , because p{x,t) < xi for (x,f) G A 2 . Moreover, since > 0 in B{ei) x I and p = Q 
in Ai U A 3 by definition, we get that ip > 0 on Ai U A 2 U A 3 . Further, by (12.31) and (I2.10L 
ipt — Alp — cip > 0 in U. Thus the parabolic maximum principle implies that ip > 0 in U. Now, 
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remember that s{x*) = + ei). By direct calculation ^{x*, 1) = 0 and, since = 0 

on dB{ei) x I, we have that 




(xM) 


1 

\/2 dei 


(xM) 


diy 


(x*,l)) = 0. 


This implies that ^(x*, 1) = 0 = ^ !)• Since ■0 > 0 in U, it follows that ^^(x*, 1) > 0, and 

therefore 


^2^61 

ds'^ 


M)>m0(xM) 


—Afie 


-0r^-t 


l-^(x -y)-s = 2ne 


which yields the first inequality in (12.91) with e := 2fxe > 0. For the second inequality, note 

that the function is antisymmetric in x with respect to Hiei), and therefore 


^2 ei 

^(-M) = ^(-M) = 

where s(x*) = + ei)- The proof is finished. □ 

The next lemma will be helpful to guarantee positivity near corner points. 

Lemma 2.5. Let w € (7^(-B(ei)) be a function such that the following holds. 

(i) w = 0 on dB(ei). 

(ii) There is e > 0 such that 


and |^<-£ foraUaeH(e,}naB. 

os[xj‘‘ ds{x)^ 

where s(x) = + si), s(x) = + ei), and v is the inwards unit normal 

vector field on dB. 

(Hi) There is a function x ■ [0,diam(i?)] —> [0,oo) with lim x('i?) = 0 such that \D^w{x) — 

■I?—fO 

D‘^w{y)\ < x(|x - y\) for all x,y £ B{ei). 

Then there exists 5 > 0 depending only one, y, and B such that w >0 in [H[ei)C\dB]sC\B(ei). 

Proof. By hypothesis (ii) and (Hi) there is (5i = 5i(£, y) > 0 such that, for every x* € H{ei)r\dB 
we have that 


^ d§'(x*)2 < -g forallxe B5,(x*)nS(ei). (2.11) 

If > 0 —where Ai is as in m — we also assume that 

We show that the claim holds for any 6 £ (0, ^]. Without loss of generality we may assume 
for the rest of the proof that the domain B is two-dimensional, since otherwise we can repeat the 


8 

















following argument to w restricted to -B(ei) fl P where P is any plane containing Rei to yield 
the result. 

By contradiction, assume there are x* G Il(ei) n dB and y G B{ei) D Bs{x*) such that 
w{y) < 0. Then there are yi G H{ei) fl Bs^ (x*) a nd Ai > 0 such that y = yi + Ais(a;*). Let 
Li := {A > 0 : yi + As(a;*) G B{ei) (1 Bsi{x*)} and A := {yi + As(x*) : A G Li}. Note 
that A n OB ^ 0, by our choice of S and (12.121) . In particular, we may find A_b > 0 such that 
2/1 + Abs(x*) G dB. Moreover, due to the fact that B is assumed to be two-dimensional, there 
is A 2 > 0 such that x* + \ 2 s{x*) =: 2/2 G A. Let L 2 := {A > 0 : x* + \s{x*) G B(ei) (~1 Bsj^(x*)} 
and define the functions fi : Li —>■ R and /2 : L 2 —t R by 


/i(A) := w(yi + As(x*)), / 2 (A) := w(x* + As(x*)). 


By (I2.11|) . we have that > e in £ 2 - By Assumption (i) it follows that 


m 


dw(x*) 

ds{x*) 


1 / dw{x*) 
y/2\diy{x*) 


dw{x*)\ 
dei ) 


= 0 = w{x*) = / 2 ( 0 ). 


Therefore / 2 (A) > 0 for A £ L 2 n( 0 , 00 ). Since 2/2 G A, there is A 3 > 0 such that fiiA^) = w{y 2 ) = 
/ 2 (A 2 ) > 0. But then /i(0) = 0, /i(Ai) = w{y) < 0, fiiAa) > 0, and /i(Ab) = 0. This contradicts 
the fact that /{' < —e in Li by (I2.11|) . Therefore w > 0 in Bs{x*) n B(ei) as claimed. □ 

We are ready to prove a perturbation result for supersolutions of scalar equations. 

Lemma 2.6. Let v G C^’^(B x /) be a function such that (resp. satisfies (12.31) and 

such that assumptions and (Hk) hold. Then there exists p > 0 depending only on a, /3o, 

k, and B such that u® (•, 1) > 0 (resp. —u® (•, 1) > 0^ in B{e') for all e' G with |ei — e'| < p. 

Proof. Assume that satisfies (12.31) and that {Ha,po) (iLfc) hold. By Lemma and 
assumption {Ha^pg), there are ei > 0 and pi > 0 such that 


r(a:, l)>ei and (x, 1) < —£1 for x G dB n H{e') 


ds(x)^ 


ds{xy 


and for e' gS" ^ with je' — ei| < pi. By Lemma 1^31 there is di £ (0, such that 


'(•,!)> 0 in ldBnH(e%,nB(e') 


(2.13) 


for all e' £ ^ with |e' — ei| < pi. Further, by Lemma [2.31 with i5 = (5i there is £2 > 0 with 

1) > £2 on dB(ei)\[dBnB(ei)]s^. Then, by {Ha,po), there are <52 £ (0,(5i) and p 2 G (0,pi) 


du 

with the property ^^{-, 1 ) > £2 in [dB{e')\[dBr\H{e')]si]s 2 for all e' £ with |e' —ei| < p 2 . 

Since u® (•, 1) = 0 on dB{e') we have 

v^' (•, 1) > 0 in [dB{e')\[dB n H{e')]s,]s 2 n B{e') (2.14) 

for all e' G with |e' — ei| < p 2 . Moreover, by Lemma 12.31 with S = 62 there is p > 0 with 

1) > pxi > ^ for all x £ B{ei)\[dB{ei)]s 2 - Then, by (iLa,/3o); there is ps £ ( 0 ,^ 2 ) such 

that 


(x, 1) > for all X G B{e')\[dB{e')]s 2 


(2.15) 
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for all e' G with |e' — ei| < ps. Therefore (I2.13p . (12.141) . and (12.151) imply that 1) > 0 
in B{e') for all e' G with |e' — ei| < p^. It is easy to check that p^ depends only on a, /3o, 

k, and B so this yields the claim for . 

From the proof of Lemmas 12.31 and 12.41 it is clear that their claims also hold if we write — 
instead of . Therefore we can argue exactly the same if we now assume that satisfies 
(12.31) . This ends the proof. □ 

We now use Lemma l2.6l to prove a perturbation result for systems. First we fix some notation. 
For a pair of functions (tii, V 2 ) G C^’^{B x /) D C{B x /) x C^'^{B x /) D C{B x I) and e G 
define the functions f: i? —>■ K. by 


vl{x,t) := Vi{x,t) - Vi{ae{x),t), n^Xjf) 

:= V2{cre{x),t) - V2{x,t). 

(2.16) 

Assume that 

- ciiu®^ - 012^2^ > 0 

in B{ei) x I, 


{v 2 '-)t - An2" - C2in®^ - 022 ^ 2 ^ > 0 

in B{ei) x I, 

(2.17) 

=vl^ =Q 

on dB{ei) x I, 

> 0 and > 0 

in B{ei) x I, 



where 


C 12 > 0 and C 21 >0 in B{ei) x I. (2-18) 

To simplify the notation, we write 

Q :=B(ei) X [i,y] and Qs := {B{ei)\[dB{ei)] 5 ) x [^ , ^], (2.19) 

for ^ > 0, where was defined in (12.21) . Note that Qs CC Q. 

Proposition 2.7. Let 7 G (0,1) and M > 0 be given constants and let vi,V 2 G C^’^{B x /) D 
C{B X I) be functions such that 

\vi\ 2 +T,Bxi < M fori = 1,2. ( 2 . 20 ) 


_ _ 2 

Suppose (ni\n 2 ^) satisfies (12.171) . (12.181) . and \\cij\\L'^{B{ei)xi) < M. 

i,j=0 

(i) If there is ^ > 0 and ei > 0 such that 

> £1 infci 2 >£i, (2-21) 

Qs 

then there exists cr > 0 depending only on ei, 5, M, 7 , and B such that ||n®^H^oo^g) > cr, 
where Q and Qs are as in (12.191) . 

(a) If there is £2 > 0 such that ||w®^ ||l°“(Q) > £2 for i = 1,2, then there is p > 0 depending only 
on £ 2 , M, 7 , and B such that 

Vi (•, 1) > 0 in B{e') for i = 1, 2 and e' G with |ei — e'\ < p. 
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Proof. To prove (i), note that (12.211) and (12.201) imply that there is a nonempty open subset 
H C Qs such that > ei in 17 and |17| > cri for some cri depending only on Si, S, 7 , M, and 
B. Let (p € Cf°{Q) be a function depending only on S and B such that ip > 0 in Q and ip = 1 in 
Qs- Then 

A := ci 2 V^^ (pd{x,t) > / ci 2 V 2 ^ d{x,t) > efai, 

Jq Jq. 

by (12.181) . (I2.21|) . and the choice of </?. On the other hand, 

j - ciivl^]ip d{x,t) = [ {-ipt - Aip - ciiip)vl^ d{x,t) 

JQ JQ 

< lkriU“(Q) / d{x,t). 

JQ 

by (12.171) and integration by parts. Therefore ||ioo(Q) > Ap. =; o- > Q, where C := Jq \p}t\ + 
|A(^| + M\p\ d{x,t). Since cr depends only ei, 5, 7 , M, and B this proves the first claim {i). For 
the second claim (ii), note that 

- Ciivl^ > Cijvf >0 in B{ei) x /, for i,j = 1,2 i ^ j, 

by (12.171) and (I2.18|) . Then, for i = 1, 2 we have that Vi satisfies the assumptions of Lemma 12.61 
with a = k = e, and /3o = M. The application of Lemma gives a constant p > 0 such that 
Vi (•, 1) > 0 in B{e') for * = 1, 2 and e' G with |ei — e'| < p, as claimed. □ 

2.1 Homogeneous linear equations 

We quote some estimates from [9] . This will be used in the proof of Theorem 11.21 to guarantee 
condition ( 12 . 211 ) for the linearized system. 

Lemma 2.8 (Particular cases of Lemma 3.9 and Corollary 3.10 in [^). Let v G C^’^(B x (0,oo))n 
C{B X [0,oo)) be a positive solution of 

Vt — Av — cv = 0 m B X (0, 00 ), u = 0 on dB x (0, 00 ), 

where ||c||loo(bx(o,oo)) < Po fov some /3o > 0. Then, there are positive eonstants Ci, C 2 , and 1 ) 
depending only on B and j3o such that 

II -> C'i(dist(a;,5B))’’, x G B, r G [1, 00 ), (2.22) 

and 

lk(-,T + t)||L=o(B) 

lk(GT)|U=.(S) 

where dist(x,9B) := inf{|x — y\ : y G dB}. 

By combining (I2.22|) and (12.231) we get the following. 

Corollary 2.9. Assume the hypothesis of Lemma \2.8l For any k > 0 there are positive constants 
C and d depending only on B, /3o, and k such that, for any r > 2k, 

C > — > C'“^(dist(x, dB))^, X G B, t G [—k, k]. 

\\v[-,T)\\l’=^(^B) 


C2 


1 


T G ( 1 , 00 ), t G [0,1], 


(2.23) 
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Proof. Let v : B x {0, oo) —>■ M be given by v(x, t) := v{x, kf) for x € t > 0. Clearly v satisfies 
the assumptions of Lemma 1^151 Let Ci, C 2 , and d be the constants (depending only on B, Pq, 
and k) given by Lemma for v. Note that (|2.2fi|) implies that 


IR-,r + t)|U=o(B) ^ C2J’ 


T e ( 1 , 00 ), t G [ 0 , 1 ], 


Then, by STTBi . (IT^ . and (EH]), 

Cf^\\v{-,T)\\La.(^B) > v{x,t + T) > C'i(dist(x,9B))’’||fi(-,t + T)||ioo(B) 
> C'iC2(dist(x,5B))’’||fi(-,r)||L=o(B) 


(2.24) 


for all X G i3, r > 2, and t G (—1,1). The result follows. 


□ 


3 Regularity of solutions 

Let a G (0,1], fc a nonnegative integer, a = k + a, and il C be a domain. Put Q := fix (t, T) 
for 0 < r < T. Following [131 page 4] we define := {/ : |/|a;Q < 00 }, where 

[/]a;Q :=sup/ : (x,t),( 2 /,s) G Q, {x,t) ^ {y,s)\, 

lF-2/r + F-sh J ,2^, 

\f\a-,Q-= sup|i:)fL>^/|+ [D^Dlf]a.Q, 

\f3\+2j<k '3 |/3|+2i=fc 


and D^,Dl denote spatial and time derivatives of order fj G and j G No respectively. 

The following is the Dirichlet version of [M Lemma 2.1] complemented with standard regu¬ 
larity arguments (see for example [T31 Remark 48.3 and Remark 47.4 (iii)] or El Theorem 2.2]) 
to guarantee (7^+'’'’^+'’'/^—regularity if the right hand side of the equation is Holder continuous. 

Lemma 3.1. Let K > 0 be a given constant, LI C be a smooth bounded domain, / C K open, 
g G L°°(Ll X /), and let v G x /) fl C(fl x I) be a classical solution of 

Vt — Ax = g{x, t) in Lt X I, v = 0 on dfl x /, 


where l|xj|Loo(Qx/) + ||sl|L“(nx/) < K. Let I C I with dist(1,9/) > 2. Then there are constants 
C > 0 and 7 G (0,1), depending only on LI and K such that 


l'^ll-|-7;nx[s,s-|-2] — ^ 


for all s £ I. 


Moreover, if lglQ;nx/ < K for some a > 0, then there are constants C > 0 and 7 G (0,1), 
depending only on LI, K, and a such that 


l'^l2+7;nx[s+l.s-|-2] — ^ S £l. 

In particular, if I = (0,oo), then the semiorbit {v{-,t) ■. t > 1} is relatively compact in C{Ll). 
Therefore uj(v) is a nonempty connected compact subset of C{Ll) satisfying lim inf j|x(-,t) — 

t->-oo z£ij{v) 

^llL“(n) = 0 . 


12 















Proof. Fix s S I and set Q := x [s, s + 2] C fl x I. Then v G because v G C^’^(Q). 

Then [TUI Theorem 7.30, p.181] implies the existence of Ci > 0 depending only on il and K such 
that 

\\D^v\\lN + 3(^Q) + ||Ut||LN + 3(Q) < C'l(||g||iN + 3(Q) + 11 U11 iiV + 3(Q) ) < 2Cl\Q\'^K. 

Then, by a standard interpolation argument, ||u||^,J,, 2 ,l < C 2 for some constant C 2 {^,K) = 

C 2 > 0. By Sobolev embeddings (see, for example, [TU| embedding (1.2)] and the references 
therein), we then have that v G for 7 = 1 — G (0,1), and there is a 

constant C 3 (fl) = C 3 > 0 such that 

|^^|i+ 7 ;Q < - C'3C'2- (3.2) 

Now, assume | 5 |a;nx/ < for some a > 0. Let 7 = min{a, 7 } and define w : B x I —>■ M 
by w{x,t) := {t — s)v{x,t). It follows that w = 0 on dpQ —here dp denotes the parabolic 
boundary — and Wt{x,t) — Aw{x,t) = {t — s)g{t,x) — v{x,t) in Q. Then, by [121 Theorem 48.2] 
(see also [TUI Theorem 4.28] and [H Section 1.6.6]), there is a constant C 4 > 0 depending only 
on n and a such that 


|wl2+7;Q < C'4(l|wl|icx,(Q) + \tg - v\^.,q). 

Then, by (]3.2L there is some constant C 5 > 0 depending only on K, fl, and a such that 
jwj2+7;Q < C5. This implies that kl2+7-r2x[s+i s+2] — ^ fo^ some constant C > 0 depending 
only on K, 17, and a. Finally, the last claim follows from the Arzela-Ascoli Theorem (see [TUI 
Proposition 53.3]). □ 


4 Dirichlet parabolic systems 

For this section let 


ui,U2 G X (0,00)) n C{B X [0,00)) n L°°{B X (0,00)) (4.1) 

be nonnegative functions such that u = {ui,U2) is a classical solution of (11.11) and such that 
assumptions (hi) and (h2) from Theorem 11.21 are fulfilled. Note that Ui satisfies {ui)t — Aui = 
gi{x,t) in B X (0,oo) with gi{x,t) := fi{t,\x\,Ui{x,t)) — ai{\x\,t)ui{x,t)u 2 ix,t) for i = 1,2. 
Moreover, by (hl),(h2), and (14.11) there is AT > 0 such that 

||l“(Bx( o,c»)) T ||5i||L“(Bx(o,c»)) ^ A" for i = 1,2. 

Then Lemma lUUT implies that |wi|a;Bx[s.s+i] < C for all s G [1, 00), 7 = 1,2 and for some a G (1, 2) 
and C > 0. But this implies that gi is Holder continuous in i? x [l,oo) for i = 1,2. Then, by 
Lemma lUTT there are 7 G (0,1) and M > 0 such that 

|mz| 2 + 7 ;Bx[ 3 .s+i] < M for all s G [2, 00), 7 = 1,2. (4.2) 

In particular the semiorbits {Mi(-, t) : t > 1} are precompact in C{B) for 7 = 1,2, and 

lim inf l|ui(-,7) - zi1]l=o(b) + ||w2 (-,7) - ^ 2 ||l°°(b) = 0, (4.3) 

t-^CO {zi,Z2)Gui{u) 

where uj{u) is as in ILI. By a standard compactness argument, the omega limit set oj(u) of 
u = (ui,U 2 ) and the omega limit set uj{ui) are related as follows: 

oj{ui) = {zi : z = (zi,Z2) G u>(u)} for 7 = 1, 2. (4.4) 
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4.1 Linearization 


The linearization of dm is the same as in [Ml Section 5]. We recall it here for completeness. 
For e G define {u\,u%) as in (12.1611 . The same notation is used if the functions do not 

depend on time, that is, for a pair z = (zi,Z 2 ) of functions Zi : i? —>■ M, i = 1,2, we set 
zl{x) := zi{x) — zi((Te(a;)) and z^ix) := Z 2 {<Je{x)) — Z 2 {x) for all a; G i?, t> 0. 

Since u = (mi,M 2 ) solves (11.11) . for fixed e G we have 

{ul)t — Aul = chul + CI 2 U 2 in B{e) x (0, 00 ), 

{u 2 )t - Au 2 = C 2 i< + C 22 U 2 in B{e) x (0, 00 ), (4.5) 

u\ = U 2 = 0 on dB{e) x (0, 00 ), 


where G L°°{B x (0,oo)) are given for i,j = 1,2, z ^ j, a; G B, and t > 0 by cfj{x,t) := 
ai{\x\^t)ui{x,t) > 0 and cl^{x,t) := cl(x,t) — ai{\x\,t)uj(ue(x),t) with 



dvfi{t, \x\,sui{x,t) + (1 


s)Ui{(Te{x),t))ds. 


Here dyfi stands for the derivative of {t,r,v) !->■ fi{t,r,v) with respect to v. Moreover, by 
(hl),(h2), and (14.11) there is M > 1 such that 


E 


|l“(Bx(0.oo)) < AF for all e G 


jN-l 


(4.6) 


4.2 Proof of Theorem 11.21 

To prove Theorem 11.21 we use the following geometric characterization of foliated Schwarz sym¬ 
metry from [M], which is inspired by a result of Brock [3]. 

Corollary 4.1 (Corollary 2.6 of [M]). LetU C C{B) and suppose that the set 

A4 := {e G : z(a;) > z(ae(x)) for all x G B{e) and z £14} 

contains a nonempty subset Af with the following properties 

(i) Af is relatively open in ; 

(ii) for every e G dAf and z £ lA we have z < z o in B(e). Here dAf denotes the relative 
boundary of Af in 

Then there is p £ such that all elements of 14 are foliated Schwarz symmetric with respect 

to p. 

Let Af := {e £ E>^~^ : w® > 0 in B{e) x [T, 00 ) for z = 1, 2 and some T > 0}. 

Lemma 4 . 2 . The set Af is relatively open in E>^~^. 

Proof. Let e £ Af. Without loss we may assume that e = ei. Then (zz®\zz20 ^ solution of 

(14.51) . and there is T > 0 such that zz®^ > 0 in B(ei) x (T, 00 ). Let M satisfy (14.611 and (14.21) . 
Then, Proposition 12.71 fzzl applied to Vi : B x [0,1] —>■ [0, 00 ) given by Vi{x, t) := Ui{x, t — T) for 
z = 1,2 yields the existence of p > 0 such that uf{-,T -|- 1) > 0 in B(e') for e' G with 

|e' — ei| < p. Since u® also satisfies (14.5L the maximum principle for systems (see, e.g., [T3] l 
implies uf > 0 in B(e') x (T -|- 1, 00 ) and therefore e' £ Af for all e' £ with |e' — e| < p. 

Thus Af is open. □ 
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Lemma 4.3. For every e £ dM and every z £ oj(u) we have zf = Z 2 = 0 in B[e). 

Proof. Without loss we can assume that dN ^ 0. Let z = (^ 1 , 2 : 2 ) £ w(u). We will only show 
that zf = 0 in B{e) for all e £ dJ\f, since the same argument can be adjusted to show that zf = 0 
in B{e) for all e £ dJ\f. Arguing by contradiction, assume there is e £ dJ\f and k > 0 such that 

\\^2\\L°°{B{e)) > k. (4.7) 


We now use a normalization argument to construct a suitable pair of functions that satisfies 
the assumptions of Proposition l2.7l (i) and (ii). Let —>■ 00 be an increasing sequence with > 6 
and such that —>■ Zi uniformly in B for i = 1,2. Define In '■= [in ~ 3,t„ + 3] C K and 

Pn ■= \\uii',tn)\\L°°(B) for n £ N, and the functions : Bx — >• R given by Vn(x, t) := IhlEiHfor 
n £ N. By Corollary 12.91 there exists 77 > 1 and 1 ? > 0 such that 


fl ^ '^n > — dist(a;, dBF 
1 


on B X In for all n £ N. 


(4.8) 


Moreover, we have that |u„|i+^;Bx[s,s+ 2 ] < C for all s £ — 2,t„ + 2], n £ N and for some 

7 £ (0,1) and C > 0. This follows from Lemma 13.II and the fact that Vn satishes 


{Vn)t - Aun = (c - ai_U2)vn in B X In, = 0 indB X In, (4.9) 


where c{x, t) := dvfi{t, |a;|, s'«i(x, t))ds for x G B and t > 0. Note that c £ x [2, oo)) 

for some /3 > 0 by our assumptions in (hi) and (14.2|) . This implies that the right hand side of 
gH) is Holder continuous. Then, by Lemma l3.II and (03), there are 7 £ ( 0 , 1 ) and Mi > 0 such 
that 


|■U2|2+7;Bx[s,s+l] + |2+7;Bx [s.s+ 1 ] < Mi ( 4 - 10 ) 

for all s G [tn — l,tn + 1], n £ N. For e £ and n £ N we also consider w® : H x /„ —>• R 

given by Vn{x,t) := Vn(x,t) — Vn{o'e{x),t), and we note that 

«)t - = Ci 2 ul in B{e) x In, 

{ul)t - Am 2 - C 22 ii 2 = C 21 < in B{e) x In, (4.11) 

= ^2 = 0 on dB{e) x In, 

where £12 '■= ^ = Q^iVn, C 21 ■= c^i/Sn = C 2 i\\ui{-,tn)\\L°°(^B), and the coefficients cf^ are as in 
(14.51) for i,j = 1,2. Then, by (14.61) . (14.81) . (hi), and (h2), there is M > Mi such that 

l|Cii||L“(B(e)x/„) + l|C22lU“(B(e)x/„) , . 

+ l|Ci2llL“(B(e)x/„) + l|c2lllL“(B(e)x/„) < M. 

for all e £ By (14.71) . (14.3|) . (14.101) . and Remark [2Tl there is (5 > 0 such that, passing to a 

subsequence of tn, we have 

\\u 2 {-,tn)\\L--{B{e)\[dB{e)]s > k for all n £ N. (4.13) 

Moreover, by assumption (h2) and (14.8|) . there is ei £ {0,k) such that 

Ci 2 = aiVn > — dist{x,dB)^ > ei in {B{e)\[dB{e)]s) x In (4-14) 
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for all e € ^ and n G N with 6 as in (03). For the choices of M, S, 7 , and £i made in this 

proof, let (T > 0 be the constant given by Proposition 12.71 (i). Without loss we may assume that 
a < k. Moreover, let p > 0 be the constant given by Proposition 12.71 (ii) with £2 = (t and M, <5, 
7 , as above. 

By (14.101) . (I4.13|) . and since e G dM^ there is e G Af with |e — e| < p/2 and 

l|u 2 (->in)|U-(B(g)\[aB(g)], > fc for all n G N. (4.15) 

Because e G Af and —>■ 00 , there is h G N such that > 0 and u? > 0 in B{e) x Taking 

into account (|4.1()L (14.1111 . and (I4.12L we see that the functions wi, ^2 : B x [0,1] —>■ 00 given by 
wi{x, t) := Vfi{x, tn — j+t) and W 2 {x, t) := U 2 {x, tn — j+t), satisfy the assumptions of Proposition 
12.71 Then, by (14.141) and (I4.15|l . Proposition 12.71 (i) yields that > <t. 

This together with (14.151) implies, by Proposition 12.71 (ii). that {-jtn + |) > 0 and {-jtn + 

|) = > 0, for all e' G with |e' — e| < p. Then the maximum principle for 

parabolic systems (see, e.g., [lU Ch. 3 Sec. 8 Thm. 13]) implies that U 2 > 0 and uf > 0 in 
B{e') X {tn + |,oo), and therefore e' G M for all e' G with |e' — e| < p. But this yields in 

particular that e G A/”, which contradicts the fact that e G dJ\f ^ since J\f is open by Lemma 14.21 

Therefore zf = 0 in B{e) as claimed and this ends the proof. □ 

Proof of Theorem \ 1. Define U := uj{ui) U —uj{u 2 ) = {zi,—Z 2 : z G a;(u)}, where the last 
equality is a consequence of (14.41) . Also let M := {e G > 0 in B{e) for all z G U}. 

Then we have that Af C A4. Moreover, for e G as in (h3), we have m®(-, 0) >0, u®(-, 0) ^ 0 
in B(e) for i = 1,2. The maximum principle then implies that u® > 0 on B(e) x (0, 00 ) for 
i = 1, 2, so that e G N and thus N is nonempty. Moreover A/" is a relatively open subset of 
by Lemma [4.21 and, by Lemma [4.31 z = z o for all z G W and e G dJ\f. The result now follows 
from Corollary 14.11 □ 

Remark 4.4. The upper bound in (14.81) holds only in the intervals = [tn — 3,t„ + 3] for 
n G N. In particular, the norm ||un(‘, Sn)||L°°(B) could tend to infinity for some sequence s„ —>■ 00 
(this would occur, e.g., if iti(-, t„) —>• 0 but mi(-, s„) 74 0 as n —>■ 00 ). In this case the coefficient 
C 12 in the linearization (|4.11l) can not be uniformly bounded independently of n G N. This is 
the reason why the eventual positivity of the perturbed difference functions—(u® , ) in the 

proof of Lemma 14.31 — is crucial. An approach based on maximum principles for small domains as 
in [5] would have the advantage of requiring weaker regularity assumptions on the nonlinearity, 
but this technique typically only guarantees asymptotic positivity of (u® , U 2 ) and therefore, 
to perform a perturbation argument, one needs that the coefficients of the linearization are 
uniformly bounded independently of (u® , ) for times, which, as explained above, does not 

happen in general. For scalar equations, eventual positivity using maximum principles for small 
domains was obtained in [S] Proposition 6.11], but an extension of this result to the case of linear 
systems remains, up to our knowledge, an open question. 
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